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Abstract The deconfinement phase transition from hadronic matter to quark matter in the 
/■^ , interior of compact stars is investigated. The hadronic phase is described in the framework 



of relativistic mean-field (RMF) theory, when also the scalar-isovector (5-meson effective 
field is taken into account. The MIT bag model for describing a quark phase is used. The 
i-^ ■ changes of the parameters of phase transition caused by the presence of i5-meson field are 

explored. Finally, alterations in the integral and structural parameters of hybrid stars due 
Q i to both deconfinement phase transition and inclusion of (5-meson field are discussed. 

t/^ . Key words: Equation of state: mean-field: neutron stars: quarks: deconfinement phase 

^1 ' transition 

> ■ 1 INTRODUCTION 

(^ , The structure of compact stars functionally depends on the equation of state (EOS) of matter in a suffi- 

lO ' ciently wide range of densities - from 7.9 g/cm^ (the endpoint of thermonuclear burning) to one order 

'^ \ of magnitude higher than nuclear saturation density. Therefore, the study of properties and composi- 

[■~~. ■ tion of the matter constituents at extremely high density region is of a great interest in both nuclear 

^^ I and neutron star physics. The relativistic mean-field (RMF) theory (IWaleckalll974l Serot & Waleckal 

^D . 1986, 19 971) has been e ffectively applied to describe the struct ure of finite nuc l ei JLala zissis et al. 19971; 

Typel&Wolter II 19971) . the features of heavv- ion collisions (IKo & Li Ill996t IPrassa et al. Il2007h . and 
the equation of state (EOS) of nuclear matter JMiller & Serot |[l995l) . Inclusion of the scalar-isovector 



J^ ■ (5-meson in this theoretical sc heme and investigation of it s influence o n low d e nsity asymme t ric nuc lear 

rS \ matter was realized in Refs. iKubis & Kutschera I ( Il999h ; iLiuetal. I ( l2002h ; lOreco et al.1 ( l2003h . At 

j^ ■ sufficiently high density, different exotic degrees of freedom, such as pion and kaon condensates, also 

deconfined quarks, may appear in the strongly interacting matter. The modern concept of hadron-quark 
phase transition is based on the feature of that transiti on, that is th e presence of two conserved quantities 
in this transition: baryon number and electric charge dOlendenning 1992) . It is known that, depending 
on the value of s urface tension, a,, the phase transition of nuclear matter in to quark matter can occur 
in two scenarios (iHeiselberg et al. lll993HHeiselberg & Hjorth-Jensen II 19991) : ordinary first order phase 
transition with a density jump (Maxwell construction), or formation of a mixed hadron-quark matte r 
with a continuous variation of pressure and density (Glendenning construction) ( lGlendenning||l992h . 
Uncertainty of the surface tension valu es does not allow to determine the phase transition scenario, tak- 
ing place in reality. In our recent paper ( lAlaverdvanl2009ah in the assumption that the transition to quark 
matter is a usual first-order phase transition, described by Maxwell construction, we have shown that the 
presence of the (5-meson field leads to the decrease of transition pressure Po, of baryon number densities 
npf and uq. 
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In this article we investigate the hadron-quark phase transition of neutron star matter, when the 
transition proceeds through a mixed phase. The calculations results of the mixed phase structure 
(Glendenning construction) are compared with the results of usual first-order phase transition (Maxwell 
construction). Also influence of (5-meson field on phase transition characteristics is discussed. Finally, 
using the EOS obtained, we calculate the integral and structural characteristics of neutron stars with 
quark degrees of freedom. 

2 NEUTRON STAR MATTER EQUATION OF STATE 

2.1 Nuclear Matter 

In this section we consider the EOS of matter in the region of nuclear and supranuclear density (n > 0.1 
fm^-^). For the lower density region, corr esponding to the ou ter and inner crust of the star, we have used 
the EOS of Baym-Bethe-Pethick (BBP) JBavm et al. Ill97lb . To describe the hadronic phase we use the 
relativistic nonlinear Lagrangian density of many-particle system consisted of nucleons, p, n, electrons 
and isoscalar-scalar (a), isoscalar-vector (w), isovector-scalar (S), and isovector-vector (p) - exchanged 
mesonslJ 



^i(a^a(x)a''a(x) - m„a{xf) - U{a{x)) + ^^mluj^' {x)u: ^{x) - U 



+4'^{iY''df,~me)ipe, (1) 

where x = x^ = {t,x,y,z), a{x), io^{x), S{x), and ~f^^{x) are the fields of the a, ui, S, and p 
exchange mesons, respectively, U{(t) is the nonlinear part of the potential of the cr-field, given by 
iBoguta & Bodmeriil977, 

Uia)^^mN{g,af + ^{g,a)\ (2) 

TOjv, irie, rric, m^, mg, nip are the masses of the free particles, ip^f = I ^^ ] is the isospin doublet 

for nucleonic bispinors, and ~t^ are the isospin 2x2 Pauli matrices. The symbol "^►" denote vectors in 
isotopic spin space. This Lagrangian also includes antisymmetric tensors of the vector fields w^ {x) and 

"^^(x) given by 

^fj.i' {x) = df_,uj^ {x) - dyuj^ {x) , SR^t, {x) = 0^"^j, (x) ~ d„~^ ^ (x) . (3) 

In the RMF theory, the meson fields a {x), lo^ {x), 5 (x) and 7^ (x) are replaced by the effective 

mean-fields W,uJf^ ,d and ^ . 

This Lagrangian density (1) contains the meson-nucleon coupling constants, g^r, g^i, gp and gs, 
also parameters of cr-field self-interacting terms, b and c. In our calculations we take for the S cou- 
pling constant, ag = (gs/mg)'^ = 2.5 fm^, as in Refs. iLiu et al. I ( |2002|) : iGreco et alTI ( l2003h : 
[Alaverdyan (2009a), for the bare nucleon mass, rrij^ = 938.93 MeV, for the nucleon effective mass, 
771^ = 0.78 771 jv, for the baryon number density at saturation, uq =- 0.153 fm^'^, for the binding energy 
per baryon, /o = — 16.3 MeV, for the incompressibility modulus, K ~ 300 MeV, and for the asymmetry 

/rjX Q O O 

energy, Eiym = 32.5 MeV. Five other con stants Oq- = (q^/mj j) , a^^ = {gui/rriui) , flp — {gp/nip) , b 
and c then can be numerically determined (lAlaverdvanl2009al) . 



' We use the natural system of units with h - 



Influence of Scalar-Isovector (5-Meson Field on Quark Phase Structure in Neutron Stars 
Table 1 Model Parameters with and without S -Meson Field 

Oct ( fm^) a^ (frn^) ag (fm^) Up ( fm^) 6 ( fm^^) c 



RMFaojpS 9.154 
RMFatop 9.154 



4.828 


2.5 


13.621 


1.654 ■ 10"^ 


1.319-10 


4.828 





4.794 


1.654- 10-2 


1.319- 10 



In Table 1 we list the values of the model parameters with and without the isovector-scalar S meson 
interaction channel (The models RMFaujpS and RMFaujp, respectively). 

The knowledge of the model parameters makes it possible to solve the set of four equations in a 
self-consistent way and to determine the re-denoted mean-fields, a = g^a, w = (?c^cJo, S = gs5^^\ and 
p = gpPo^^\ depending on baryon number density n and asymmetry parameter a = (n„ — np)/n. The 
standard QHD procedure allows to obtain expressions for energy density e{n, a) and pressure P{n, a) 
of nuclear npe plasma: 



A;_ {n,a) 



eNM{n,a,pe) ^ — / Jk^ + m*{a,S)'^ k'^dk 



1 







jWTm*JaJ^ k^dk 


-- 771JV fT^ + - (T-* + - — + — + — + i— 

6 4 2. \aa a^ as ap 



1 



v/fc2 + me2 k^dk, (4) 





k- {n.a) 



Pi\!M{n,a,pe) = — / (Jk-{n,ay+m*{(T,S)'^-Jk^+77i*{a,Sy^jk'^dk 


^ ^ ^^/k+in, a)2 + m* (a, S)^ - ^k^ + m*„{a, 6)A k^dk 



^2 







b 3 C 4 1 / ^2 C^2 ^2 2 

6 4 2 \ Qa ttui as Up 



\/Mf-mf 



-—^Pei^J-l-ml) I \/W+m^ k'^dk, 

OTT^ TT J 



(5) 





where fie is the chemical potential of electrons, 

TO*(o', S) ~ m.N — a — 6, m^{a, S) ~ niN — cr + S (6) 

are the effective masses of the proton and neutron, respectively, and 

Q 2 \ 1/3 



k±in,a)=[^—il±a)\ . (7) 
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The chemical potentials of the proton and neutron are given by 

^p(ri, a) = Jk^ {n, a)"^ + m*{a, 5)"^ +oj + -p, 



Unin, a) = y/k+{n, a^ + m;^{a, S) ^ +uj - -p. 



(8) 
(9) 
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Fig. 1 The asymmetry parameter as a func- 
tion of the baryon number density n for a 
P -equilibrium charge-neutral npe -plasma. 
The solid and dashed curve correspond to the 
RMFoujpS and RMFaojp models, respec- 
tively. 
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Fig. 2 Re-denoted meson mean-fields as a func- 
tion of the baryon number density n in case 
of a /3 -equilibrium charge-neutral npe-plasma 
with and without 5-meson field.The solid and 
dashed curves correspond to the RMFcjujp5 
and RMFaojp models, respectively. 



In Fig.l we show the asymmetry parar neter a for the /3 -eq uilibrium charge-neutral npe-plasma, as 
a function of the baryon number density, n d Alaverdvanl2009al) . The solid and dashed curve correspond 
to the RMFauipS and RMFaujp models, respectively. One can see that asymmetry parameter falls 
off monotonically with the increase of baryon number density n. For a fixed baryon number density n, 
the inclusion of the J-meson effective field reduces the asymmetry parameter a. The presence of (5-field 
reduces the neutron density n„ and increases the proton density np. 

In Fig. 2 we plotted the effective mean-fields of exchanged mesons, a, u, p and (5 as a function of the 
baryon number density n for the charge-neutral /3-equilibrium npe-plasma. The solid and dashed lines 
correspond to the RMFaiopd and RMFaojp models, respectively. 

From Fig. 1 and Fig. 2 one can see that the inclusion of the scalar-isovector virtual 5{aQ (980)) meson 
results in significant changes of species baryon number densities Up and ?i„, as well as the p and 5 meson 
effective fields. This can result in changes of deconfinement phase transition parameters and, thus, alter 
the structural characteristics of neutron stars. 

The results of our analysis show that the scalar - isovector (5-meson field inclusion leads to the 
increase of the EOS stiffness of nuclear matter due to the splitting of proton and neut ron effective 
masses, and also the increase of asymmetry energy (for details see Ref. lAlaverdvanll2009al) . 

2.2 Quark Matter 



To describe the quark phase an improved version of the MIT bag model JChodos et al. 119741) is used, in 
which the interactions between u, d and s quarks inside the bag are taken into account in the one-gluon 
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exchange approximation dFarhi & Jaffe Ill984l) . The quark phase consists of three quark flavors u, d, s 
and electrons, which are in equilibrium with respect to weak interactions. We choose to„ = 5 MeV, 
nid = 7 MeV and iris = 150 MeV for quark masses, and as = 0.5 for the strong interaction constant. 

2,3 Deconflnement Phase Transition Parameters 

There are two independent conserved charges in hadron-quark phase transition; baryonic charge and 
electric charge. The constituents chemical potentials of the npe-plasma in /3-equilibrium are expressed 
through two potentials, yuj^ and ^^^ , according to conserved charges, as follows 

M«=a4"^^'' Mp = /^i^^^^ - /^ef ^^\ /ie=Mef^^^- (10) 

In this case, the pressure Pnm, energy density enm and baryon number density tinm, are functions of 
potentials, /i^ ' and /ij.j ' . 

The particle species chemical potentials for udse-plasma in /3-equilibrium are expressed through 
the chemical potentials /ij, ' and ^]J ' as follows 

1 / (QM) _n {QM)\ 



Md - Ms - o (^A^fc +Me/ j: (11) 

(QM) 
Me = Mrf - Mu = Mef ■ 

In this case, the thermodynamic characteristics, pressure Pqm, energy density eqm and baryon 

number density riQ A/, are functions of chemical potentials /i^ and ^^ . 

The mechanical and chemical equilibrium conditions (Gibbs conditions) for mixed phase are 

M^^^Mr^'^M., Mir=Mir^=Me. (12) 

Pquit^b, IJ'el) = PNMiy'b, Me/)- (13) 

The volume fraction of quark phase is 

X = Vqm/ [Vqm + Vnm) , (14) 

where Vqa/ and Vnm are volumes occupied by quark matter and nucleonic matter, respectively. 

We applied the global electrical neutra lity condition for mixed quark-nucleonic matter, according 
to Glendenning JGlendenning II 1 992l l2000l) . 



2 11 

- UuiHb, MeO - o "rf(Mfc, Me;) - o "s(Mfo: MeO - "-e(Me/) 



(1 - X) ["p(Mf), MeO - neiHel)] 

0. (15) 



The baryon number density in the mixed phase is determined as 

"- = (1 -X) [npi^l'b, ^lel) +nn[^J.b,^J.el)] 



- xlnuip-b.pi-ei) + nd{^ib,^J■el) + nsinb, ^J■el)] , (16) 



and the energy density is 



£ = (1 - X) [£pi^J'b, MeO + e«(M&: Mez)] 

+X l^uifJ'b, Me;) + ed(M&> Me;) + es(M&> Me;)] + £e(Me;)- (17) 
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Table 2 The Mixed Phase Threshold Parameters with and without 6 -Meson 
Field for Bag Parameter Values, B = 60 MeV/fm^ and B = 100 MeV/fm^ 



Model 



njv 
(fm-3) 



riQ 
dm-'') 



(MeV/fm-^) 



Pq 

(MeV/fm^) 



(MeV/fm^) 



(MeV/fm^) 



B60(Ta;p(5 


0.0771 


1.083 


0.434 


327.745 


72.793 


1280.884 


B60(Ta;p 


0.0717 


1.083 


0.336 


327.747 


67.728 


1280.889 


BlOOcrcjp^ 


0.2409 


1.448 


16.911 


474.368 


235.029 


1889.336 


BlOOo-cjp 


0.2596 


1.436 


18.025 


471.310 


253.814 


1870.769 



In case of x = 0, the chemical potentials jjl^ and /i^, corresponding to the lower threshold of a 
mixed phase, are determined by solving Eqs. (13) and (15). This allows to find the lower boundary pa- 
rameters Pjq, Ejy and njy. Similarly, we calculate the upper boundary values of mixed phase parameters, 
Pq, £q and rig, for x = 1- The system of Eqs. (13), (15), (16) and (17) makes it possible to determine 
EOS of mixed phase between this critical states. 

Note, that in the case of an ordinary first-order phase transition both nuclear and quark matter are 
assumed to be separately electrically neutral, and at some pressure Po, corresponding to the coexistence 
of the two phases, their baryon chemical potentials are equal, i.e.. 



i^NM {Pq) = iJ.QM {Po) 



(18) 



Such phase transition scenario is known as phase transition with constant pressure (Maxwell construc- 
tion). 

Table 2 represents the parameter sets of the mixed phase with and without (5-meson field. It is shown 
that the presence of (5-field alters the threshold characteristics of the mixed phase. For B = 60 MeV/fm^ 
the lower threshold parameters, njv, sn and Pn, are increased, meanwhile the upper ones, nq, eq and 
Pq, are slowly decreased. For B = 100 MeV/fm'' this behavior changes to opposite. 
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Fig. 3 EOS of neutron star matter with the de- 
confinement phase transition for a bag constant 
B = 60 MeV/fm . For comparison we plot 
both the Glendenning and Maxwell construc- 
tions. Open circles represent the mixed phase 
boundaries. 



Fig. 4 As in Fig 3, but for B = 100 MeV/fm^ 



In Fig. 3 and Fig. 4 we plot the EOS of compact star matter with deconfinement phase transition 
for two values of bag constant, B = 60 MeV/fm^ and B = 100 MeV/fm'^, respectively. The dotted 
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curves correspond to pure nucleonic and quark matters without any phase transition, while the solid 
lines correspond to two alternative phase transition scenarios. Open circles show the boundary points of 
the mixed phase. 
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Fig. 5 Constituents number density versus 
baryon number density n for B — 100 MeV/fm 
in case of Glendenning construction. Vertical 
dotted lines represent the mixed phase bound- 
aries. The dashed curves show appropriate re- 
sults of the model without J-meson field. 



Fig. 6 As in Fig 5, but for Maxwell construc- 
tion. Vertical dotted lines represent the density 
jump boundaries. 



In Fig. 5 we plot the particle species number densities as a function of baryon density n for 
Glendenning construction. Quarks appear at the critical density tim = 0.241 fm^'^. The hadronic mat- 
ter completely disappears at ng — 1.448 fm^'^, where the pure quark phase occurs. The solid curves 
correspond to the case, when also the 5- meson effective field is taken into account besides a, lo, p 
meson fields (model i3100_ aujpS). The dashed curves represent the results in case when we neglect the 
(5-meson field (model i3100_ aup). One can see that inclusion of the (5-meson field leads to the increase 
of number densities of quarks and protons, and simultaneously to the reduction of number densities of 
neutrons and electrons. In Table 2 we have already shown the mixed phase boundaries changes caused 
by the inclusion of the 5 - meson effective field. 

Fig. 6 shows the constituents number density as a function of baryon number density n for B = 100 
MeV/fm-^, when phase transition is described according to Maxwell construction. Maxwell construction 
leads to the appearance of a discontinuity. In this case, the charge neutral nucleonic matter at baryon 
density ni = 0.475 fm^'^ coexists with the charge neutral quark matter at baryon density n2 = 0.650 
fm^-^. Thus, the density range ni < n < n^ is forbidden. In case of Maxwell construction, the chemical 
potential of electrons, pe, has a jump at the coexistence pressure Pq. Notice, that such discontinuity 
behavior takes place only in usual first-order phase transition, i.e., in the Maxwell construction case. 

3 PROPERTIES OF HYBRID STARS 

Using the EOS obtained in previous section, we calculate the integral and structural characteristics of 
neutron stars with quark degrees of freedom. 

The hydrostatic equilibrium properties of spherical symmetric and isotropic c ompact stars in 
general relativity is described by the Tolman-Oppenheimer- Volkoff(TOV) equations JTolman I Il939l : 
lOppenheimer & Volkoff |[l939l) : 



dP _ G {P + e){m + 4:Trr^P) 
~ih ^ "^ 1 - 2G m/r ' 



(19) 
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dm 
~d^ 



= Airr £, 



(20) 



where G is the gravitational constant, r is the distance from the center of star, m{r) is the mass inside 
a sphere of radius r, P{r) and e{r) are the pressure and energy density at the radius r, respectively. 
To integrate the TOV equations, it is necessary to know the EOS of neutron star matter in a form 
e(P). Using the neutron star matter EOS, obtained in previous section, we have integrated the Tolman- 
Oppenheimer Volkoff equations and obtained the gravitational mass M and the radius R of compact 
stars (with and without quark degrees of freedom) for the different values of central pressure, P^. 
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Fig. 7 The ma.ss-radiu.s relation of neutron 
star with different deconfinement phase tran- 
sition scenarios for the bag constant B — 
60 MeV/fm^. Open circles and squares denote 
the critical configurations for Glendenning and 
Maxwellian type transitions, respectively. Solid 
circles and squares denote hybrid stars with min- 
imal and maximal masses, respectively. 



Fig. 8 As in Fig 7, but for B = 100 
MeV/fm . The mass-radius relation in case 
of the Glendennig construction without the 5- 
meson effective field is also displayed for com- 
parison (dash-dotted curve). 



Fig. 7 and Fig. 8 illustrate the M{R) dependence of neutron stars for the two values of bag constant 
B = 60 MeV/fm'^ and B = 100 MeV/fm'^, respectively. We can see, that the behavior of mass-radius 
dependence significantly differs for the two types of phase transitions. Fig. 7 shows, that for B = 60 
MeV/fm"^ there are unstable regions, where dM/dPc < between two stable branches of compact stars, 
corresponding to configurations with and without quark matter. In this case, there is a nonzero minimum 
value of the quark phase core radius. Accretion of matter on a critical neutron star configuration will 
then result in a catastrophic rearrangement of the star, forming a star with a quark matter core. The range 
of mass values for stars, containing the mixed phase, is [0.085Mq; 1.853Mq] for B ~ 60 MeV/fm'^, 
and is [O.QOTA/©; l.TSOA/©] for B = 100 MeV/fm^. In case of MaxwelHan type phase transition, the 
analogous range is [O.2I6A/0; I.828M0] for B = 60 MeV/fm''. From Fig. 8 one can observe, that in 
case of i? = lOOMeV/fm^, the star configurations with deconfined quark matter are unstable. Thus, 
the stable neutron star maximum mass is 1.894M0. Our analysis show, that for B = lOOMeV/fm'^, 
the pressure upper threshold value for mixed phase is larger than the pressure, corresponding to the 
maximum mass configuration. Hence, in this case, the mixed phase can exist in the center of compact 
stars, but no pure quark matter can exist. The dash-dotted curve in Fig. 8 represents the results in case 
when we neglect the (5-meson field (model i?100_ autp). One can see that for a fixed gravitational 
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Table 3 Hybrid Star Critical Configuration Properties for i? = lOOMeV/fm^ 
with and without (5-Meson Field 





Minimum Mass Config 


:uration 


Maximum Mass 


Configuration 


Model 


(MeV/fm^) 


Mmin 

(Mo) 


R 
(km) 


(MeV/fm^) 


Mmax R 

(Mq) (km) 


BlQOaujpS 
B100aa;p 


235.029 
253.814 


0.997 
0.955 


14.354 
13.960 


1390.77 
1386.03 


1.780 11.190 
1.791 11.139 



mass the star with (5-meson field has larger radius than the corresponding star without the (5-meson field. 
Influence of (5-meson field on the hybrid star properties is demonstrated in Table 3, where we display the 
hybrid star properties with and without (5-meson field for minimum and maximum mass configurations. 
The results show that the minimum mass of hybrid stars and corresponding radius are increased with the 
inclusion of flie ^-meson field. Notice that influence of (5-meson field on maximum mass configuration 
properties is insignificant. 

4 CONCLUSIONS 

In this paper we have studied the deconfinement phase transition of neutron star matter, when the nuclear 
matter is described in the RMF theory with (5-meson effective field. We show that the inclusion of 
scalar isovector (5-meson field terms leads to the stiff nuclear matter EOS. In a nucleonic star both 
the gravitational mass and corresponding radius of the maximum mass stable configuration increases 
with the inclusion of the 5 field. The presence of scalar isovector (5-meson field alters the threshold 
characteristics of the mixed phase. For i? = 60 MeV/fm'^, the lower threshold parameters, n^r, £jv; Pn, 
are increased, meanwhile the upper ones, nq, eq, Pq, are slowly decreased. For B = 100 MeV/fm^ 
this behavior changes to opposite. 

In case of the bag constant value B = lOOMeV/fm^, the pressure upper threshold value for mixed 
phase is larger, than the pressure, corresponding to the maximum mass configuration. This means that 
in this case, the stable compact star can possess a mixed phase core, but the density range does not allow 
to possess a pure strange quark matter core. 

Stars with (5-meson field have larger radius than stars of the same gravitational mass without the 
(5-meson field. Alterations of the maximum mass configuration parameters caused by the inclusion of 
(5-meson field is insignificant. 

For the bag constant value B = 60 MeV/fm'^, the maximum mass configuration has a gravitational 
mass Mmax = 1-853 Mq with radius R = 10.71 km, and central density pc = 2.322 • 10^^ g/cm'^. 
This star has a pure strange quark matter core with radius tq sa 0.83 km, next it has a nucleon-quark 
mixed phase layer with a thickness of rMP ~ 9.43 km, followed by a normal nuclear matter layer with 
a thickness of r^ « 0.45 km. 
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